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AERONAUTICAL SYMBOLS. 
I. rCNOAMENTAL AND DiiRIVEU UNITS. 





Sycubal. 


Length- 


I 


Time.... 


t 


Force 


F 


Power. . . 


V 


Speed — 





Metric !. 



Unit. 



mfitor 

MiM'ond 

vvoi^rlit of one kilogram. 



mh,(H^ 



Symbol. 



111. 
see. 
kg- 



j ' 

i Englisli. 

1 


Unit. 


Symbol. 


! foot (or mile) 


ft. (or mi.). 


' ^erond (or hour) 


sec. (or hr ). ! 


wc;ight of one pound 


i 




IP 

M. V, H, 



2. GENERAL SYMBOLS, ET(\ 

Weight, ^V= m<h Specific wei^lit of ^ '.standarcr* air, 1.223 b^m.' 

Staiidard aocM^lr»raUvfn of gravity, — ().070o5 ll>/rt.^ 

^ = 9.S06m/.se^.''^--C2.L72Lt/;i.ec.^ Moment of inertia, rnlc'^ (indicate axis of the 

ir racJiiis of ^^yration, t, by proper subscript), 

f/ Arci), *S'; winj^ area. iSV, etc. 

Density (ma^s per unit volume), p Gap, G 

Standard density (if dry air, 0.1247 (kg.-m.- vSpan, h; chord length, c, 

sec.) at 15.GT. aad 700 mm. -0.00237 (lb.- Aspect Yii\v^^llc 

ft.-sec.) ])istance ivom c, g. to elevator hinge,/. 

OoefHcient of visco.sity, ^. 



Mass, m- 



True airspeed, V 

]\vnamic (or impact) pre^i^aro, ^7 = ,!^ pT 

L 



3. AERODYNAMICAL SYMBOLS. 

Dihedral angh^, 7 



Reynolds Number = p— , where Z is a linear di- 



Lift, L; absolute coenicient 
Drag, D; absolute coofhrient 



D 



Cross-wind force, C: ar):^aiute co{^(Ueient 
C 



qS' 



nicnsioii. 

e. g., for a niodrl airfoil 3 in. chord, 100 mi/hr., 
nornnd pressure, OX': 255,000 and at 15.G°C, 

230,000; 

or for a model of 10 cm. chord, 40 m/.sec., 
corresponding numbers are 299,000 and 

270,000. 

Center of pressure coefficient (ratio of distance 
(Note that these coeflicients are twice as of C. P. f mm leading edge to chord length), 
large as the old coefilcients Lcr ^K-^' (-r- 
Angle of setting of wings (relative to thrust Angk> of -tabilizcr setting with reference to 

line), iw lower wing. 

Angle of stabilizer setting with reference to Angle of at i nek, a 
thrust line it Angle of downwash, c 



Resultant force, R 



t 
! 
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MATHEMATICAL EQUATIONS FOR HEAT CONDUCTION IN THE FINS 

OF AIR-COOLED ENGINES. 

By D. R. Harper 3d, and W. B. Brown. 



SUMMARY. 



At the request and with the support of the Engineering Division of the Air Service United 
States Army. McCook Field, the Bureau of Standards undertook some laboratory nVe^^^^^^ 

Committee on Power Plants for Aircraft and by that committee recommended for puWicat on 
as a technical report of the National Advisory Committee for Aeronautics. In connec on 
with laboratory measurements of tk. heat-dissipating power of typical engine cylinder a 
mathematical analysis of fin behavior was made and is given in this commu5cation 

Ihe introduction contains a description of the paper which wUl assist the general reader 
who 18 not mterested in mathematical detail in finding those parts of the paper most likelv 
to prove useful to him. A recapitulation of the mathematical developments is given in 
bection IV and forms the statement of conclusions reached in so far as a mathematical paper 
of this type may be said to have conclusions. Numerical examples illustrative of these con- 
equltfons''^^ * ''^'^ ^"'^ suggestion of possible appUcation of the 

The problem considered is that of reducing actual geometrical area of fin-cooling surface 
which IS of course not uniform in temperature, to equivalent "cooling" area at one definite 
temperature, namely, that prevailing on the cylinder wall at the point of attachment of the 

»-«n ^^A^T ^ '•'^'^"g ^'^rf*^ ^ if it P-^t of the cylinder 

wall and 100 per cent effective. j 

The quantities involved in the equations are the geometrical dhnensions of the fin thermal 
conductmty of the material composing it, and the coefficient of surface heat dissipation between 
the fin and the air stream. Several assumptions of physical nature are thus necessarUy involved 
'IToh Z^ the problem possible of solution. These are set forth in detail, and the lunitations 
pointed out ^PP^y^^S the equations to numerical calculations are carefully 

An expression for approximate fin effectiveness is developed, based upon simple mathe- 
matics and very convenient in form for engineering use. The essence of the paper is an 
examination into the magnitude of the errors involved in using this expression without cor- 
rection and a determmation of the corrections needed for accurate work. The mathematical 
expressions mvolved are quite complicated, including Fourier's Series, super Fourier's Series 
Bessel functions of zero order of two kinds with imaginary arguments, etc. The results of 
the work are collected in graphical form in a series of charts, so that the design engineer can 
use the simple formula first developed and apply to it corrections readily read from the 
charts, thus avoiding entirely all higher mathematics. 

I. INTRODUCTION. 

The equations which express the flow of heat in a metal in terms of simple physical proper- 
ties are perfectly definite and adapted to numerical computations, although usually somewhat 
cumbersome and tedious. In applying these equations to the fins on the outside of the cylin- 
ders of au--cooled mternal-combustion engines, the chief obstacle to numerical work is the 
great uncertainty of the value to be assigned to one important pliysical quantity, the rate of 
dissipation of heat from the fin surface into air under the conditions surrounding the fin. This 
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important deficiency seems to have discouraged any widespread use of the equations of heat 
conduction in considering the problem, since deductions made from them could be trusted 
only within rather wide limits. , v 

With the increasing knowledge of rate of cooling m an air stream, it has become more 
worth while to compile the information obtainable from a mathematical analysis of the prob- 
lem The details of such an analysis are not of sufficient general interest to warrant the average 
reader in foUowmg them closely, but it is believed desirable to render them available for refer- 
ence by those who are working in the same field. The equations which have to do with this 
subject are bulky, the algebra and integrations tedious and tune consuming, and the chances 
for error are high, although no especially intricate or abstruse reasoning is involved, nor is there 
much difficulty in interpreting final results other than the necessity of a careful exammation of 
the relative numerical size of the various terms. - , , ■ t t 

With a full appreciation of the tax unposed on a reader by reason of the foregoing facts, 
the authors have prepared this paper in a form which endeavors to meet the following speci- 

fication: , i i_- ^ r ■ 

(a) To segregate, for the benefit of aU who are L^iterested in the general subject of air- 
cooled engines, a general skeleton of the analysis, includmg ihe discussion of conditions which 
bear upon the problem, statement of the exact assumptions to which are applied the mathe- 
matical development and the conclusions resulting therefrom, with a few examples of numerical 
computations to illustrate the practical application of the mathematics. 

(6) To interline with the above, in such form that the general reader may skip it without 
losing the thread of the development, such details of mathematical transformation as will be 
needed by the specialist to reproduce the equations or use them to advantage m theu- applica- 
tion to his particular problem. 

(c) To omit all detaUs of algebra, integration, arithmetic, etc., which are merely the 
mechanism of the mathematical development. Although these steps are essential to an 
acceptance of the validity of any of the deductions, it must be left to the critical reader to supply 
the gaps, because the paper is sufficiently complicated in meeting specification (6) without 

Biich additional weighting. . ■ ^i, • j- 

The basic principle of design which characterizes an air-cooled engme is the providmg 
of some means to increase greatlv the natural surface by additional cooling surface, the purpose 
being to keep the engine cylinder wall temperature dowTi to a value below the upper limit 
set by satisfactory engine performance. This additional surface takes the form of cooling 
fins usually made an integral part of the cylinder barrel and arranged either longitudinally 
alon«- the barrel or circumferentially around it. The problem considered is that of reducing 
actutl geometrical area of cooling surface, which is, of course, nonuniform m temperature, to 
equivalent "cooling" area at a definite, easily specified temperature. This may be done by 
finding an expression for the effectiveness of fin surface, i. e., the ratio of the amount of heat 
dissipated by unit area of fin surface to that dissipated by an equal area of cylmder wall surface 
with the same tempearture as that at the fin base. This will make it possible to treat aU the 
cooling surface as if it were on the cylinder wall and had 100 per cent effectiveness. 

II. APPROXIMATE FORMULA FOR EFFECTIVENESS OF COOUNG FINS. 

Two cases occur in practice: (1) Circumferential fins, usuaUy of tapering thickness, with a 
base temperature that may change from point to point; (2) longitudinal fins with sunilar con- 
ditions of thickness and temperature. A direct analytical investigation of each case m all its 
generality is quite difficult and has not yet been completely worked out. In this paper an 
indirect attack on the general problem is made by a method of successive approxunations. 
The effectiveness is first computed for a simple case, where several simplifying assumptions 
are made. Then the limitations of these assumptions are removed, one or two at a tune, and 
the necessary correction made to the first result. The effectiveness / is, therefore, expressed 
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GENERAL ASSUMPTIONS. 

Four general assumptions of physical nature are made that apply to every case considered 
here : 

(1) Quantity of heat transferred per unit time from the metal surface to the air is pro- 
portional to the temperature difference between the metal and the air, 

(2) The coefficient g, heat transferred in unit time from a unit surface per unit temperature 
difference, is constant over the fin surface. 

(3) The fin is symmetrical about a plane through its middle and approximately parallel 
to its faces. 

(4) The temperature at a given point is independent of the time. 

Assumption (1) is known as Newton's Law of Cooling. It has been found to be sensibly 
true for very high velocity air-stream cooling in such cases as are under consideration.* 

Some preliminary measurements of heat dissipation by air-cooled engine cylinders made in 
the laboratories of the Bureau of Standards also indicate the validity of the assumption within 
limits necessary to work of this kind. 

Assumption (3) is obviously true mechanically except for minor inequalities of manufacture, 
which on an average over a cylinder would be inappreciable. Physically, there is marked 
lack of symmetry about such a plane if the fin be oblique to the air stream, but this lack of 
symmetry has to do with assumption (2). 

Assumption (4) provides a working basis that is entirely acceptable. Questions of engine 
design must be settled from considerations of cooling capacity under full load and continuous 
operation, when a steady temperature distribution exists. A cooling system which will meet 
this demand will obviously meet the less stringent cooling requirements involved in starting 
up and approaching temperature equilibrium. The rapid variations in temperature of the 
inside wall of the cylinder between explosion and intake are quite damped out in their effects 
on fin temperature. The validity of the assumption has been demonstrated by experimental 
evidence.* 

Assumption (2) involves, amongst other things, either independence of the heat transferred 
from a fin to an air stream and the velocity of the air stream, or else the assumption of constant 
velocity of the air stream over all portions of the fin. The first hypothesis is untenable; the 
sec<md one is discussed below. Experimental data are very incomplete. The assumption is 
recognized as weak, but in the present state of knowledge it is about all that can be done. It 
is known to be justifiable in long tubes and is probably not far wrong for longitudinal fins ^ith 
fairly open spacing. One set of measurements ■ on a plain cylinder without fins, a small cylin- 
drical rod (diameter 2 cm., i inch) indicated that with the air stream normal to the axis of the 
cylinder, the variation in air velocity, front and rear, was of the order of 30 per cent, namely, 
a difference of 15 per cent each way from the mean value. British measurements on air velocity 
between fins indicate a change in q from tip to root of longitudinal fins of less than 15 per cent. 
These meager data would suggest that if an average value over a fin were taken, the deviations 
from it would not generally be more than 10 or 15 per cent. The error after integration in such 
cases is generally less than the original error, so that the result with this approximation is 
probably in the neighborhood of the true value. 

SIMPLIFYING ASSUMPTIONS. 

In addition, the following simplifying assumptions give a simple problem which serves as 
a basis for more complete analysis: 

(1) The temperature across the fin thickness is sensibly constant. 

(2) The fin is so long that the effect of the exposed ends (of longitudinal fins) is neglible. 

» A. H. Gibson. Automotive Industries, May 13, 1920, p. 1109. Theories and Practices in the Air Cooling of Eugiaef, 
" Judge, "High-Speed Internal Combustion Engines," pp. 107-109. 

• T. E. Stanton, Great Britain Advisory CommlttM for Aeronautics, Report 94, 1912-13, p. 47. 
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(3) The heat loss from the exposed edge can be accounted for by imagining the width 
extended by a distance equal to one-half the fin thickness at the outer edge and assuming no 
heat loss from the end. 

(4) The base temperature is constant, 

(5) The fin thickness is constant. 

(6) The fin is longitudinal. 

The details of assumption (3) are illustrated in Figure 1, 



a H 

34 



The edge abc at a temperature 6, dissipates some heat A^j^. If ah and he be swung around 
to ah' and ch' and the space filled in with metal, the temperature drop from a to ¥ will be small 
and these surfaces will dissipate practically the same amount of heat as before. The other 
simplifying assumptions require no explanation. 

It will be shown later that all of the A fs are small compared with so that the cross 
products which represent the effects due to the joint action of several disturbing influences 
may be neglected and only first-order effects due to single causes need be worked out. 

OUTLINE OF THE GENERAL PLAN OF MATHEMATICAL ANALYSIS. 

First,/' is computed on the basis of all six simplifying assumptions stated above. 

Second,/ is computed with assumptions (1) and (3) removed. Thus/-/' gives A^/^f 
^tZ+^ijs/ These are shown to be negligible for all conditions prevaiUng in engine work. 

Third, /is computed with assumption (4) removed, giving A^/ which is zero. 

Fourth, / is computed with assumption (5) removed and replaced by the one that the 
sides are straight, and hence the thickness t is a linear function of the distance x from the fin tip. 
This gives Ag/, which is often not negligible. If \ is the ratio of the thickness at the tip to the 
average thickness, then Ag/ may be expressed in terms of/' and X. Its value is computed for 
several values of X and the results shown graphically on a chart. 

Fifth, / is computed with assumption (6) removed and a longitudinal fin replaced by a 
circumferential one. This determines A«/ If p is the ratio of the inner radius to the outer, 
/ may be expressed in terms of/' and p only, so that the results for A^/ are expressed graphically 
in the same way as above. 

By means of these charts, / can be found with a sufficient degree of approximation for the 
most general case. 



^ ^ 1 










\ 
t 
1 




» J^' \ 





(1) Computation of/'. (The approximate function to represent heat flow in terms of 
geometrical and thermal properties of the fin.) 

Let 5 = coefficient of heat transfer from surface. 
w = true width of fin. 
^ = fin thickness. 

x = distance from the cylinder wall. 

/=fin effectiveness. 

t^fin conductivity. 

^ = temperature of fin above the air. 
^0 = temperature of cylinder wall above the air. 
i7 = heat dissipated by the fin per unit time, 
fffl heat dissipated by equal area of wall surface per unit time. 

w' ^ corrected fin width t<;+ 
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By introducing the simplifying assumptions (3) and (4), the boundary conditions are: 

^=*0 when x«=i0' 
whena:«0 

'rhe fundamental equation of heat transfer* under the conditions corresponding to assump- 
tions (1), (2), and (5) is & F 

d^e 2qe 

d^^ld (2) 
The solution of this equation in the form most convenient for this work is 

B^A cosh a (x~B) (3j 

where A and B are arbitrary constants of integration and a is an abbreviation for V2?/R 
^^When A and B are determined to satisfy the boundary conditions, stated mathematically 

^ = ^ <^osh a (x — w^) 
® cosh anv' 

The rate of heat dissipation from a unit length of fin is computed by an integration with 
respect to x from 0 to w\ 

B-2j\edx^ J;' cosh a (X - v>^) & (4) 
The heat dissipation from an equal area {2w') of cylinder wall at temperature $^ is 



and therefore 



This function (tanh aw')/aw\ for which the single letter / will be used, is the function 
which will serve as the basis of discussion for much of the following work. Under average con- 
ditions of practice it will be found sufficiently exact to serve as the basis for all computations. 
In those cases where it does not fit with sufficient accuracy, it will be found convenient to use 
it as a principal term plus necessary correctmg terms. The principal properties of the function 
therefore merit attention. 

When aw' is increased, the value of tanh aw' increases likewise, but rather slowly, and 
although reaching 0.9 when aii)' = 1.50, it does not increase beyond 1.0, no matter how large 
aw' becomes. The ratio therefore, starts at unity and gradually decreases to zero, when 
plotted against values of aw\ This plot is shown in Figure 3, where the single letter u is sub- 
stituted for aw'. 



(6) 



< The method of deriving this equation is explained fully in elementary textbooks on heat. The difference in quantity of heat conducted into 
an elemental slab at coordinate r and that conducted out at z+dt la, in the equfUbritlili state, the amount that escapes into the air tbroueh the 
surface dx. This equality, when common factors are canceled, is tbe equation (2). 
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Figure 3 shows that / grows less as u increases, whence as or q increase, or as or < 
decrease. That is to say, the effectiveness of a fin decreases as — 

(1) The fin width is made greater. 

(2) The fin is made thinner. 

(3) The fin is made of material with a poorer thermal conductivity. 

(4) The coefficient of heat dissipation to the air increases. 

This last statement indicates that fin effectiveness is relatively less at higher air speeds, 
since the value of q increases with air speed. 



/oo 



^60 



^60 

I 



-ix /son abbreviafion for aw w 'h-^ 
w'-rfn widih-i- */z fin thickness, 
- k*Therm<^con<kjctMfyof fin, 

_q ^ Coefftcienf of surface 

heaf dissipafion, 
t - Thickness of fin. 




Fig. 3.— Curve showing function 



\/a/u€S of 12 



2.0 a.^ 



2.8 



Any units whatever if they 



belong to a mutually consistent system will do to measure the above 
four quantities and w^ill lead to the same number for f 

III. CORRECTION TERMS TO APPROXIMATE FORMULA. 

In the use of the formula (5) derived in the preceding section, it is very necessary to know 
the departure from exactness which has been introduced by use of the simplifying assumptions 
which are obvious deviations, at least to a small degree, from the conditions which actually 
prevail. This investigation of the order of magnitude of the several correction terms involves 
considerable tedious mathematics which are incorporated in the paper only to meet the needs 
of a specialist in the field. The general reader will find the conclusions summarized at the end 
of the chapter. 

(I) CORRECTIONS FOR HEAT DISSIPATION FROM EDGE OF FIN AND FOR "CROSS-FLOW" IN FIN. 

This correction is determined by solving the problem as stated initially without the aid of 
sunplifying assumptions (1) and (3) but with the aid of the remaining simplifications. The 
difference between the solution so obtained and the function defined in equation (5) above 
shows the value of the correction which would be necessarj^ to the latter to take account of the 
two factors of this title. 



Heaf , 



FigA 



The problem is a straightforward development of two-dimensional heat flow', in which the 
boundary conditions are stated mathematically as follows: 

With the origin of coordinates located as indicated in Figure 4, the axis y = 0 is the median 
line of the fin, and from considerations of symmetry there can be no heat flow across the median 
plane, a condition expressed mathematically by zero temperature gradient. 



» Byerly-Fourlers Series and Spherical Harmonics, art. 59, p. 102, edition of 1902. 
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ho 

Along 2/ = 0 g^=0 (7) 

Along y=^/2 tg^- (8) 
(this face being one which is dissipating heat to the air). 

Along x = w? ^^i'^ "9^ (9) 
Along ^ = ^0 (10) 

(this edge of the fin being the one integral with the cylinder wall and maintained at given con- 
stant temperature ^o)- The fundamental equation of heat transfer for this problem is 

which, being of the second order in two variables, can include in its primitive four arbitrary 
constants and permit of applying the four conditions (7) to (10) for the determination of auch 
constants to give a complete solution of the problem as stated. 
The convenient form for the primitive is 

oo 

O^^XmAm cosh a (x-B) cos ay (12) 

which contains three arbitrary constants A, B, a and satisfies the condition (7). The develop- 
ment of condition (8) leads to determination of a as any solution of 

ctt . at at .^^ 
2 tan2=|fc (13) 

Then the value of B is defined in terms of a by the use of condition (9) leading to 

a tanh a (5 - ti;) =« I (14) 

Lastly, values of Aja must be selected to conform to condition (10) in terms of the a and jB 
already specified. These values of A must satisfy 

BQ^XmAm cosh aB COS ay (15) 

The possibility of determining values of A^^ to satisfy this relation has been established by 
workers in the field of Fouriers Series and other harmonic expansions,' and while the particular 
form here appiied may not identify exactly with those commonly found in the textbooks, it 
seems quite unnecessary to supply here any of the intervening transformations. It is adequate 
evidence of the validity of the assumption that definite values for -4 m may be found if we pro- 
ceed to find some which define a convergent series. 

The principal reductions will have to do with equation (13) and the algebra is simplified 
considerably by the use of a parameter 4> in place of a, defined by the relation 

at 

(16) 



• Bytrly-Fouiisn Series and Spherical Harmoiilcs, pp. U8~121. 
3341—23 2 
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It may be noted that equation (13) is one with an infinite number of real values of a or * 
satisfying it, whereas equation (14) has but one value of B for any given value of a, the hyper- 
boHc tangent being a single valued function. To dispose of this equation and replace it by a 
series, the quantity B is replaced by a new parameter €, defined by 

B = + « (17) 

In terms of t and * the revised equations (13) and (14) take the form 

«tan«=|j (18) 

With a determined as a root of equation (18) and the corresponding €m determined from 
the series (19), we must select an to give the expansion (15) which in the new parameters is 

^0= Sm^« cosh ^(w -f 2 + cos (20) 

and then substitute all these in (12) revised to 

(?=Sm^„ cosh cos (21) 

The details pertaining to the process just outlined are tedious and of interest only to the 
worker who desires to check the development. Multiply each side of (20) by cos (2«ky/<) dy, 
where *k is any root of equation (18) and integrate for y from 0 to i/2, 

J\ cos dy=sm^^ cosh T^(^+1+*™)JJ' ^-^y ^y (22) 

nr^ 26m 2(6k 7 tr^m sin </>m cos </>ic-</)k sin <^»k cos 

eos-^l/cos-p-2/(/2/=-2L --<^-<t>^'~^ J 

and therefore vanishes for those values of 0m and 0k which are roots of equation (18) except 
for the particular choice <t>k = <t>my where the indeterminate form of the expression introduces 
the possibility of fmite value. It follows, therefore, that the summation on the right-hand 
side of equation (22) can consist of no more than the smgle term given by the equality of m 
and Jcj whence 

r«,c>»?f! i, dn-A^ o«h ^-f (w + J" ^t" y iy 
=^„co*[^|=(„.|...)]x|[i.5£*.] 



From which 



A^-\ sech ^ + 1 + *«) J J cos ^ y dy (23) 

sin 20m 



1+- 



20„ 



The above value of is in a form for a general expansion of as any arbitrary func- 
tion of y, namely, for any specified temperature along the thickness edge of the fin (i. e., along 
the cylinder wall). This generality is too complicated for consideration in this paper, which 
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has been limited to the case of 6^ being constant. This condition may therefore be intro- 
duced at this point into equation (23), which reduces to 

A^^4e, sech ?fB («, + ,/2 + .„) (24) 

This value of used in equation (21), with the aid of equations (18) and (19), completes 
the solution of the problem as a problem in thermal conduction. 

For application to the purposes of this paper, it is desired to compute the effectiveness of 
the fin surface, which has been defined above as the ratio of the heat dissipated by the fin to 
the heat which would be dissipated in the same time by an equal area of surface all maintained 
at the constant, uniform temperature 6^, which is here the temperature of one edge of the 
fin. Designate this ratio as H/Hq, and retain the previous notation with q for coefficient of 
surface heat dissipation, remembering that 6 is measured on a temperature scale with its zero 
chosen at the temperature of the air stream cooling the fin. 

H^^q (2w + 0 »o H^qjedS 

To integrate over the dissipating surface, we have dS = dx along the two fin surfaces for 
unit length of fin and dS = dy along the edge opposite the cylinder wall. 



jedS = 2 J\-t/2 dx + 2 J J ^x- w dy 



I dx -h ^x-w dy 

H^+l) 

From equation (21) 

CO / f \ 

tfy«t/2 = Sm^m cosh ~ (x-W-^- Cm J COS 0„ 

<^x-w =2m cosh 2±a (^ + e„) COS y 
and by substituting these values in equation (25) and reducing the result, 

20„ 



. t ^ sm2</»„j 

J 4^m 



^_^t I'^iP^ (2«m + sm2*J 
2 



tanh 



t °° 



o «:«3 JL cosh 
2 sin' </>, -* 



*m (2.^m+sin 2«„) 2^„ ^ , t , -\ (26) 



In the use of equation (26) the terms of the summation are given by using for <j>ai in succes- 
sion all the real roots of equation (18) 

0 tan = (18) 

g = heat dissipated from fin to air stream per unit time per unit area of fin surface per 

unit temperature diflference. 
^=«fin thickness. 

^ = thermal conductivity of fin material. 
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From each successive value of <t>m, the corresponding is computed from equation (19) : 



It is obvious that equation (26) is much too cumbersome for everyday application. From 
the standpoint of pure mathematics it is the equation to choose in preference to equation (5) , 
being based upon more acceptablo assumptions and including as a special case the conditions 
leadhig to the solution expressed by equation (5). From the mathematical viewpoint, correct 
procedure would be to solve the problem in the more general form, obtaining equation (26) 
for effectiveness of a fin, and then to examine the relative magnitude of the terms involved 
and leave out the small ones in order as successively less exact approximations are desired. 
In this way one might approach a comparatively simple formula, either equation (5) or an 
equivalent, for all usual applications. In this paper an alternative presentation has been 
adopted in deference to the algebraic complexity of the processes in the more general case, 
and it has been deemed wise to develop first the comparatively simple solution which is good 
enough to apply to most air-cof)lcd engine cylinder work, showing how good or how bad is 
the approximation by a later development of the more exact relations. 

These more exact relations are too complicated for convenience. The course outlined 
above can be put in practice, not by a general consideration of the magnitude of each term in the 
series, but only by specific numerical computation of these terms for one or more typical sets 

of conditions. m. t . 

The units selected for q, Jc, and t are immaterial so long as they are consistent. Ihe factor 
qt/Jc, which enters the computation, is dimensionless and independent of the unit system. If it 
be desired to measure i in inches and q in Btu. per minute per square inch per degree Fahren- 
heit, then it is only necessary to express k, the thermal conductivity, in Btu. per minute per 
square inch per unit temperature gradient in degrees Fahrenheit per linear inch. The inter- 
national units are used below. Assume: 

= 0. 10 calorics per second per cm.' per unit gradient in °C. per cm.) . 
2= 0.008 (calories per second per cm.' per °C.). 
( = 0.5 cm. 

to = 3.0 cm. .,•,/. 

This set of values is for a steel fin (fc = 0.10) of excessive proportionate thickness (one-sixth 
of the width) and large absolute dimensions, and in a very high speed wind stream, {q = 0.008 
probably corresponds to a wind velocity of 90 meters per second, 200 miles per hour.) These 
conditions are the ones which should exaggerate the effects of "cross flow" and dissipaUon 
from the thin edge of the fin, the two factors included in the more comphcated solution and 
omitted in the simpler. By "cross flow" is meant taking into account the two-dimensional 
flow of heat in a plane section of the iin rather than treating it as linear flow from the engine 
cylinder wall toward the outer edge of the fin, sensibly parallel to the flat fin surfaces and with 
but a negligible component perpendicular to these surfaces or across the fin. The case selected 
is therefore unfavorable to the approximate equation and should indicate the largest corrections 
to it necessary in any ordinary application. 

11=0,008X0.5^^^20 (27) 
2k 0.20 

The values of t^m are given by the roots of — 

</> tan 0 = ^ or tan <j> = (28) 

which is conveniently solved by graphical means by plotting values of tan <l> and of l/50<t> and 
reading the value of <p at the intersections, namely, where the two functions equal each other 
(fig. 5). The value of 1/500 are so small for all values of 4> greater than that corresponding to 
the first few roots that the function may be taken as coincident with the axis and cutting the 
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tangent curves at mw. The first root may be read from the graph of figure 5 plotted to a scale 
sufficiently open or may be approximated analytically by the series expansion, 



tan« = « + y + ||V- • 'if<t><~ 



(neglecting high powers of <i> because it is noted that 0 is much smaller than unity) , 

0.1414 = «[l+f + 'jf]* 



J5 




.05 



Fio. 5. Graphic solution of tan x 

To solve this it is evident that 4> is so closely 0.1414 that the latter may be substituted for 
0' and 4^ without appreciable error in comparison with unity. 

0. 1414 = <^ [1 + 0.00667 + 0.00005] * 
= [1 + 0.00336] 
<i> = 0.1414 X [1 - 0.00336] = 0.1409 
To find the second root, t/)2, the relation tan (ttH- a) tan a is the key to the process, bearing 
in mind that is so very close to tt as to admit of some approximations which are very exact 
while at the same time very simple. Omitting details 

02 -IT + 0.00635 

The departure of from 2ir and of 4>a from 3ir, etc., is quite inappreciable, whence the values 
of 4> for equation (26) are: 

01 - 0. 1409 02 = 3. 148 03 = 6.28 0^ 9.42, etc. 
The next step in computation is evaluation of em from equation (19). 



= 0.020 



Ar^ O.lOO 
gr"* 0.008 



= 12.50 



qt 0.020 ^-.-Q qt ^0.020_ 
2-fc0, ^ 0409 ^^^^ 2^,- 3TI48 - ^'^^^^ 



2t03 



'0.0032 



€, + 1=12.50 (0.1419)»|^l+|(.1419)' + i (0.1419)^+- » .J = 0.25356 
€,+1=12.50 (0.0063P [^1 + 1 (0.0063)'+ • • J =0.00049 
€j + |- 12.50 (0.0032)^1^1+^ (0.0032)»+- • .J = 0.00013 
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In substituting the values just obtained into the terms of equation (26) it is worth while to 
tabulate the intermediate steps for the first term (m—l) in order to illustrate the order of mag- 
nitude of the different factors and to form a guide for estimation of the changes to be expected 
in using somewhat different values for k orqoTt from those selected for this example. 



0.1409 €1 + 2 = 0.25356 

2<^>, = 0.2818 (^| + 6i) = 0.1429 

sin 2(^, 0.2780 sinh^* Q + = 0. 1434 
sin »<^i- 0.0197 cosh 



1.834 



1.0102 



, t 0.1409 (0 



0^2780 r 0.1434 -1 

^09X0.5598) ^^^^ 3.21 TO. 



1^ 4- -^ + e,- 3.2536 

tanh ^"^^ (if + 1 + €,) = 0.950 

cosh ^'(io+|+«,)= 3.21 
0.0394 1.01021 



. _ < 1^ si n 2,i>^ 

(20niH-sin20 

m/ 



t 



tanh 



1409(0.5698) 3.21 

24>, 



sinh' 



t 



(1 + *-) 



cosh 



2«„/ t , \ 



2 sin' 4f>n 



cosh 



^_^<7 4>m (2</»„^ + sin2«J 



'cosh~^^t/;+|+€«j 

/-0.1538{3.524 [0.950-0.04471 + 0.4995x0.315} +2, etc. 
«0.1638{3.189 + 0.1572}+S, etc. 

= 0.5148 + 0.1538 3,^48^6.30?^^^^^^^ ^^^^^ ^^.59 (3.0005) 
sinh 12 .59 (0.00049)"] 2 sin^ (ir + 0 .00635) cosh 12.59 (0.0004 9)) 



(29) 



49)-] , ■ ^ . 

cosh 12.59 (3.0005) J 3.148x6.31 cosh 12.59 (3.0005) j 



(30) 



Reduction of the terms in the above expression (30) indicates by inspection the general 
trend of each term in the later series. Cosh 12.59 X 3 is enormous and the two cosh terms in 
the denominators above and in all succeeding terms are so extremely large with respect to anj^ 
thing occurring in the numerators that the terms containing them are inappreciable. Also, 
the tanh function of a large argument is unity, whence there is left, to evaluate, only 

sin (2w+ 0.0127) , ~ sin 2(t>m 

3.148 (6.30 + 0.013) » 4>m (2«„+ sin 2«J 

where <t>m for m=3 and greater is so near a multiple of r that sin 2<l>n will be inappreciable. 
The numerical term above is 

0.0127 



and 



3.15x6.31-0-00064 
/= 0.5148 + 0.1538 X 0.0006 = 0.5149 



(31 
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riiis value of/ should now be contrasted with the approximate value/' given by the simple 

c>;))rcssion, equation (5): 



2=0.008 
ir^O.lO 

< = 0.50 
ii;=3.00 
to' = 3.26 



ai/>' = 1.836 
tanh aw' - 0.950 

The difference between/' and the more exact value of the effectiveness which is given by 
/ in equation (31) is therefore not quite 3 parts in 500, or less than 0.6 per cent. 

It thus appears that the error introduced by neglecting the transverse temperature gradient 
and assuming the edge correction to be simply t/2 added to w is less than 1 per cent in this exag- 
gerated instance. In ordinary cases it is negligible entirely, for an aluminum fin, g = 0.003 
(air speed 70 or 80 mi./hr.), t==0.50, u;' = 2.0 cm, / = 0.15 cm, /=/', within less than 0.1 percent. 

(2) CORRECTIONS FOR VARYING BASE TEMPERATURE AND EXPOSED ENDS. 

In this proof the edge correction that has been proved to be sufficient, i. e., w + </2 for w 
and I -\- 1/2 ^ t/2 for the length, will be assumed and the transverse temperature gradient 
neglected. Hence, the problem becomes that illustrated in Figure 6. 



•x X'ttf' 

fig. 6 

The plate of width w is replaced by the fictitious one of width w ' with the origin at the 
outer or free edge, so that the plane x^w' becomes the engine cylinder wall maintained at a 
given temperature 0^, assumed constant as to time but not uniform with relation to the coor- 
dinate 2. By neglecting the transverse (y direction) temperature gradient, the fundamental 
equation for the problem takes the form 

which must satisfy boundary conditions at the four edges as follows: 
when 

"'■ = 0 

dz' 

-0 



0 


be 




l'~ 


l+t 


0 




dx 


w' 


$ 
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where F(z) is a prescribed function, given for any particular problem. The boundary con- 
ditions for the two flat surfaces of the fin, which are cooled by the air stream, have been incor- 
porated into equation (33) in the process of deducing it from fundamental considerations. 
A convenient form of solution of equation (33) for this application is 

e-=I,A cosh OCX cos fiz (34) 

where a and P are connected by the relation 

«'-^-|? (36) 

because in the form selected for (34) the boundary conditions for x = 0 and 2 = 0 are automatically 
satisfied for all values of A and a, leaving these two parameters to be determined by the remain- 
ing two conditions. The requirement B,t z = V is satisfied if 

(36) 

where m is any integer whatever. The solution is complete in the form 

= Sm^„ cosh (^g+ 3^ COS ^ (37) 

provided values of A„ are determinable, so that when x = «?' the function (37) identifies with 
the given F(z). This is a common Fouriers Series development and requires for the value 

, - I F(z) COB ^dz (38) 



2 p'. 



Zt I J . §11, H M U " ^ 

cosh 



It is not worth while for purposes of this paper to assume any of the more probable forms for 
F{z) and complete the solution of such special cases. It happens that the general conclusion 
which is desired, namely, the difference between the approximate value of fin effectiveness which 
is given by equation (6) and the more exact value given by equation (38), is capable of being 
found in terms of a general form for Fiz) unreduced to special cases. 

To compute the effectiveness of the fin, proceed according to the definition to find expressions 
for the heat actually dissipated and that which would be dissipated if each portion of the fin were 
at the temperature of its contiguous cylinder wall. The latter quantity of heat is: 



F{z) dz (39) 
and the actual dissipation is 

dx\ edz (40) 
Jo 

From the ratio of these, the effectiveness /is 

I Odx dz 

j-sr -4^ 

"^'Jo ^^^^^^ 

When the value of 6 defined by equation (37) is substituted in (41), it contains the expressioix 

r*' mirz J 
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which is zero except for the value m = 0, and, accordingly, the working expresaion for equation 
(41) is much simplified. The zero term of a cosine Fourier expansion is half the formula for the 
general term; that is to say, equation (38) reduces to 



and 




The substitution in equation (42) of the expanded form of reduces the result to the form 




which expression is identical with the value given by equation (5), derived on the hypothesis 
that the fin base did not vary in temperature along its length. In other words, the nature of the 
variation of temperature along the base of the fin, including a uniform distribution as a special 
case, is immaterial, in so far as the fjunction to express fin effectiveness is concerned. 

(») WEDGE-SHAPED FINS. 

The complete solution for a wedge-shaped fin involves rather complicated mathematics, the 
terras involved being Bessel functions with imaginary arguments. The practical form of appU-- 
cation of the solution, expressed in manageable form for numerical work, is to plot a curve or 
series of curves giving the correction to be applied to a simple expression for the effectiveness of 
the nearest equivalent parallel-sided fin. A digest of the mathematical treatment follows. 




Let the fin be wedge-shaped, as shown in Figure 7, with straight sides, the thickness tapering 
from a value »t the fin base to a value at the fin tip. The fin width may be considered ex- 
tended in amount tJ2 to account for the heat dissipation actually occurring from the surface 
at the end (ij) and the fictitious end is then treated as if blanketed completely. The origin of 
coordinates is most conveniently located at the fictitious end and the problem is thus stated in 
terms of surface dissipation along the two sloping surfaces of Figure 7, a blanket at x = 0; namely, 
dd/dx = 0 in that plane and 6 = 0^^ at x=^v)-\'tJ2. The fin is to be considered so long (direction 
perpendicular to plane of the paper) that the end dissipation is immaterial. 

For derivation of the fundamental equation, consider the heat flow per unit time in a 
section dx at the point x, the corresponding fin thickness (0 being defined in terms of the angular 
parameter a shown in Figure 7. 

<««, + 2(a;-|) tana. (44) 



5---«2 tan a. 
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The heat flow, per unit time per unit length of fin, at the plane x is 

(where h is thermal conductivity of the fin, B the temperature at point x and x and t as defined 
in fig. 7) SO that the difference between the heat flow into an elementary slice at plane x and 
out at plane x 4- da: is 

d~x dx ) ^^"^"^^^ + ^ dx di 

By means of equation (44) this expression becomes 

U^,dx + 2ima>Jcf^dx (45) 

In the equilibrium condition, this quantity of heat equals the heat dissipated by the two 
elements of fin surface, namely, 

dx 

2qe 



cos a 



(g = heat dissipation per unit area per unit time per unit temperature difference, fin surface to 
surrounding air; temperature of fin at point x, as above, the scale of temperature used being 
such as to have zero for the temperature of the air) 

This equation is not one for which a solution may be recognized readily, the term causing 
trouble being the x in the coefficient of the second derivative. By a change of independent 
variable, involving considerable algebra, the equation reappears in a well-known form, similar 
to Fourier's equation, Bessers equation of zero order. The clue to the necessary change of 
variable is furnished by examining equation (46) in the standard form with unity for the initial 
coefficient, 

d^e 2 tana dS 2gft^ ^^^^ 

^"*"«j(l- tan a) 4-2x tan adx icosa[t, (1 - tan a)-f 2x tan a] 

do 

and trying substitutions that will shnplify the coefficient of • (It may be noted at this point 

that since in equation (46) 6 has been expressed as a function of only one variable, x, it is un- 
necessary to distinguish further between partial and total derivatives.) The change of variable 
is defined by 

where 

k sm a 

and the equation resulting from the transformation of (47) is 

p.+k^.e=0 (49) 
apt fx an 

which would be Fourier's equation were the terms all positive. Upon substituting (ifi) for/i, 
the equation goes into this form, whence the solution of (49) , as given, is 

e^AJo {ifi)^BK, (ifi) (50) 
where A and B are arbitrary constants and Jo and Ko are Bessel's functions of the two kinds, 
order zero. Before discussing briefly any properties of such functions relevant to this applica- 
tion, it is well to examine the more finished form taken by the equation when A and B are 
determined to satisfy the terminal conditions mentioned above, 

^-=0 whenx = 0 and ^=-^o when x=«;' (51) 
ax 
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These do not correspond to sunple algebraic expressions when expressed in the variable m- 



so that for all finite values of n the derivative d^/dpi vanishes whenever dB/dx vanishes, whence 

do ^ . OIL K (1-tana) 



tan . 



(52) 



In order to apply these conditions to the determination of A and B in equation (50), it is 
necessary to make use of the properties of Bessel functions which have been discovered by 
mathematicians working in the field. Even a brief review of such properties is far beyond the 
scope of this paper.' 

The Bessel function of the second kind, designated above as the Z function, is, for a com- 
plex variable, generally replaced by a slightly different form of function than that for which 
the symbol K is common in mathematical literature. The common use of K makes it the 
function related to J, so that 

(x)=logx^„ ^^^-'t-Wl)(^^'2y^(i-^'2H) 

Accordingly, with an imaginary argument (ix), a complex relation would result, 

Zo(ix)^logt/o(ix)+logxJo(ix)-^,-^2^ . . 

and it is convenient to take the term in log i over on the left-hand side of the above equation 
and define a new function which it will be noted is a real function in x. It is beyond the scope 
of this paper to consider in any detail the properties of such functions and the reasons for select- 
ing particular forms as the elements in which to express solutions. Unfortunately, there is 
considerable difficulty in comprehending the literature on the subject because great confusion 
occurs in tlie notation employed by different writers. The original extensive treatment of 
complex Bessel functions is due to HankeP and the symbol H is common for such functions 
but with exasperating lack of unanimity in regard to the exact definitions of such functions, 
which in the hands of various writers differ by several additive constants or constant multi- 
pliers. For the purposes of this paper, the H function employed will be that tabulated by 
Jahnke and p]mde,* defined by the following series: 

'/o{ix)log^^+|^ + (l+l)(2^, + (l+^ + i)^^-^ . .jwhere log 1 = 0.11593 

which it will be noted is a real function of x, since Jo (ix) is always real. 

Rewriting equation (oO) in terms of this particular form for a second solution, 

6 = AJoM + BiHoiiM) (53) 

from which 

^^^-AU,{ifji)-{-BH,M 



2 

iHo (ix) = 

IT 



dfx 



1 Among standard tqxts on Bessel flections may be cited: N. Nielsen, Handbudi der Theorie der Csrtinderfuiiktioren-Teubner, 1904- Gray 
and Matthews, Treatise on Bessel Functioitfe, MacMfllan, 1896; W. E. Byerly, l oiirier's Series and Spherical Harnjoiiics, Chop. Vll', Ginn 1902: 
Jahnke and Emde, Funktioiitafcln, Section XITI, Teubner, 1902. ' ' 

• Hankel, Mathematische Annalen 1, p. 483, l.%9; 8, p. 463, 1875. 

» Jahnke and Emde, Funktlonentafcln, p. 134 of IW edition. 
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by virtue of the properties of Bessel functions that d Jo (x)/«Zz= -J, (x) ordH^ (a)/dx= -H^ (x). 
In order that equations (52) may be satisfied, 

bJ-^^ (54) 
and A must fulfill the condition e=e„ when m=m,. 

^ „ g, MJoHm) - Ji (iMi)go(^ M) (55) 

For the particular purpose of this section, it is not necessary to tabulate numerical values 
of equation (55) and map the temperature distribution in the fin. The object sought here is 
an expression for fin efiPectiveness, defined as in the preceding sections. The heat which would 
be dissipated by the fin shown in Figure 7 to air at temperature zero if the fin surface were all 
at temperature would be (per unit length of fin per unit time) 

^ COS a ^ 

while that actually dissipated is 



dx 



cos a 

from which it follows that fin effectiveness / is 



I dx (56) 



From equation (48) 

There are two ways of integrating (57), which are in principal identical, and of course 
lead to the same result. One is to substitute for O/Oo the value given by equation (55) and 
integrate the resulting expression by using as a reduction formula 

^ [xJiix)] = xJoix), and likewise for H, 

and the other results from noting the identity 

d/ d9\_ dW de 

iiK^diir^dii^^ri^ 



From equation (49) 
whence 



m , dJB 



do 

Also, from (52), jjj vanishes for m = whence 

^"12^6, dfi]r2h'w' F,(i/i,)Jo(iM,)-e7,(tMi)Ho(W 

In computing numerical values with equation (58), the values of Mi and fi^ may be obtained 
from equations (52) in terms of w\ t,, and a, the geometrical constants of the fin. Where the 
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fin extends to a sharp-edge intersection of its two sloping surfaces, ti=^0 and therefore — 0 

and /i2"==26 V^' = 2V5'^7^ sin a. The function is infinite for zero value of the argument, 
whence for = 0 it is necessary to consider (58) in the form given below to avoid an inde- 
terminate 00 /qo. 

This special case^ a fin extending to a sharp edge, could, of course, be solved very much 
more simply than by making it a special case of the more general problem. Recurring to 
equations (49) and (50), if the condition of solution to be met is d?/dju = 0 when /i = 0, familiarity 
with Bessel functions shows at once that an abbreviated form of equation (50), namely 

e = AJo{iiJi) 

will meet the condition stated and at the same time have one arbitrary constant left so that the 
boundary condition sit x^w may be fulfilled. This method of procedure confirms very easily 
the result reached in equation (59). 

For the interpretation of equation (59), substitute for Mz i^is value 26 V^^, bearing in mind 
that for all engine cylinder fins the taper is so small that sin a = a = tan a within the accuracy 

which is required by engineering practice. The value of a is then approximately i t^/w, and 

if we denote by <m the mean thickness of the fin which tapers uniformly from at one edge 
to zero at the other, then a = t„^/w. (In this case, and w are identical.) The reduced expres- 
sion for equation (59) is 

^_ iJ. (i2.yl£) 

and expresses the eflfectiveness of a wedge-shaped fin in terms of its physical and geometrical 
characteristics. 

Following the general plan of this paper, the next step is to ascertain the difference between 
/, computed by the exact formula (60) and a value /' which would result from employing the 
ver}^ simple expression (5); in other words, assuming that a wedge-shaped fin of moderate 
taper functions nearly enough like a parallel-sided one to permit of using the formula developed 
for parallel sides and then making a correction for the error introduced by this procedure. 

tanhaw' , I^Q , ^ 

/ =-"a¥'"- ^^^^^ V 

Rewriting /in terms of a, 

in which it is seen that / is an explicit function of the product aw Just as is/', so that the two 
functions may be compared at any desired point. The diflPerence is about 6 per cent for aw='l; 
about 16 per cent for aw = 2; and 25 per cent for aw = 3. The fin dimensions commonly em- 
ployed are such that aw probably never exceeds 2 aiid is pretty generally less than 1. A com- 
parison of the two functions is plotted as Figure 8. 

It is reasonable to suppose that in all cases a trapezoid section fin would differ less in be- 
havior from that of the corresponding rectangular section fin than would a triangular section fin 
of equivalent mean thickness and width, so that in using the simple expression /^ for computing 
effectiveness of wedge-shaped fins, the maximum error occurring would be that corresponjung 
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to the difference between the two curves of Figure 8. It is therefore possible to use the simple 
formula and apply a correction with a fairly accurate degree of approximation ^ estimating the 
correction from the relative approach to a rectangle or a triangle of the fin section under con- 
sideration. A more exact but less convenient procedure is to compile a table or set of curves 
giving the exact correction under various circumstances. Such a set of curves forms Figure 9. 
The derivation is as follows : 

The taper of a fin may be expressed in terms of the fin width and the ratio of its thickness 
at the tip (<p fig. 7) to its mean thickness ^m. 

Let 




o 4 .8 is 20 e,d 

Values of aw ' 

Fig. 8. Comparison of the two functions which express effective- 
ness of a triangular section and of a rectangular section fin. J Is 
the functional for Bessel's function of the primary type, i Is the 
complex symbol V-^- Js'with an imaginary argument is a real 
function, and Ji is a pure imaginary, so that iJi Is a real 
function. 

Then a of Figure 7 may be expressed in terms of Wy t^, and X -and with no approximation 
other than sin a = tan a, substitution in the expressions (52) which define Mi and ^2 lesad to 

"-v^.H^-V^^ "-r^.^^' ' '''' 

where <i has the value used in all the previous developments, namely, ^/2q|itm' Since tm/^ is 
always small with respect to unity, 

and the term in brackets can usually be omitted. There is a 1 to 1 correspondence between 
values of aw and approximate values of fin effectiveness, so that for any desired fin effectiveness 
a suitable value of aw may be read from the curve of Figure 3, and by means of the relations 
(63), (64), Ml and M2 iiiay be tabulated as functions of X for any effectiveness. Values of and 
so obtained are then substituted in equation (58) and a comparison between the resulting / 
and the approximate effectiveness/' will give the corrections, as a function of the taper ratio X, 
which must be applied to the approximate function/'. Such curves are plotted*^ in Figure 9. 

!• The complementary procedure is to take stated values of X and for a series of such values determine the correction as a function of approximatn 
effectiveness. This procedure has been adopted for Figure 10, where it is only necessary to show the curves X«=0, X-O.,") and X=0.75 to permit of 
sufflcieatly accurate interpolation, by inspection, of any other curve of the family for the purpose of obtaining the correction to f ' for any taper at 
any effectiveness. 
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4. CIRCUMFERENTIAL FINS. 



For a circumferential fin of considerable width on the average size engine cylinder, it is not 
to be expected that relations developed for a long, rectangular fin will hold without appreciable 
correction. It will be shown here that the magnitude of the correction is well within the limits 
which justify the procedure of employnig the approximate formula, corrected when necessary, 
in preference to using an exact solution of this problem with its attendant complications. The 
difference in the physical behavior of a given area of circumferential fin and of rectangular 
plate is perhaps most easily pictured by focusing attention on the mean circumference. When 
the fin width is not small with respect to the radius of curvature of this mean, there is going to 
be a considerable difference between the fin area, for a given length of this median, which is 
within the mean circumference from that area outside it, whereas in a rectangular plate, the 
median bisects the area. With such a difference in area distribution, it is clear that the tem- 
perature of the median is not going to be by any means equal to that which is found at the 
median of the rectangular plate. The use of the rectangular plate formulas is therefore more 
in the nature of analogy than of approximation, but it is nevertheless convenient. 
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Fig. 9. —WP^lge-shapod fins. Corrections to approximate effectiveness 
Tanh aw' 



A '^^ 

Appro>itmaHs effectiveness /' 

Fig. 10— Wedge-shaped fins. C orrectlons to f. 



A further picture which may assist in visualizing the physical processes involved comes from 
comparing the way in which a circumferential fin differs from a corresponding straight one to 
the way in which a tapering fin differs from its analogue of uniform thickness. In the latter 
case, we have practically identical surfaces with differences in the metallic conducting area 
from root to tip of fin. Since the metal near the tip is less tiseful, removing a certain fraction 
there and adding it correspondingly at the root where it is most needed results in a fin somewhat 
more effective thai) the same average thickness fin with no taper. If, now, we take a straight 
fin of uniform thickness and bend it aroimd an arc, we do nothing to alter the metallic conduc- 
tion process, but do change the disposition of surface. We get a less proportion of the surface 
in near the engine cylinder, where the temperature head is greater, and a correspondingly 
greater fraction out at the rim, where it is less useful. Accordingly, the curved fin is slightly 
less effective than its straight analogue. The possibility of certain similarities in the mathe- 
matical treatment of the two corrections, taper and curvature, thus presents itself. As a matter 
of fact, it turns out that the mathematical functions involved are practically identical, although 
leading to corrections in opposite directions, as the above picture indicates. The taper cor- 
rection which has just been developed in detail and found to be always a positive correction to 
f is paralleled by one for curvature, always n^ative. 

The notation for the circumferential fins is as follows: 

fic = inner radius (namely, the outer radius of the engine cylinder). 

fif = outer radius (extreme fin radius). 
i = im thickness, assumed uniform. 

w^Rt- lie (fin width) , 

r, $, coordinates (polar) of any element of the fin. 
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The other symbols employed retain the same significance as in previous sections. The 
heat dissipated from the edge of the fin is to be taken into account by the method developed in 
section (2) of tliis paper, adding tl2 to the fin width. The exact correction, to give a surface iden* 
tically half that of the edge, would be 



but the sum of the terms following 1 is usually less than 1 per cent, and 112 is in itself only a 
small correction term, so that the omission of these terms causes no appreciable error. 

The fundamental equation of heat transfer in a metal, expressed in plane polar coordinates, 
combined with the condition for surface dissipation from both sides of each element of fin 
surface, is , ^ 

-5?-^rdr^^ ^^^^ 

{6^ temperature at any point; 5, coefficient surface heat dissipation; ky thermal conductivity; 
t, thickness, assumed uniform) . The boundary conditions are 

-|^^-Owhen r=i?'r = i?r+| (66) 

e^^Oo when r = Rc 

Equation (65) is similar to (49), in fact identical with it if fx/a be substituted for r, whence 
the form of solution is given by equation (50) or equation (53), and since the boundary condi- 
tions (66) are identically those of equations (52) with appropriate values of ar instead of Mi 
and M2» it is imnecessary to discuss any details of solution of equation (65). The .result may 
be taken by inspection from equation (55), 

H, {iaR't) Jo (iar) - (iaR't) Hp (iar) 
^-^^i?, (iaR't) Jo {iaR.)-J', HdR'f) Ho (mRc) ^ ^ 

Following the usual procedure for expressing the effectiveness of the fin surface, divide 
the heat dissipated by the fin, namely. 



rut' 



Jit' 

e'2'Kr'dr 



by that which woiild be dissipated by an equal area of cylinder wall at temperature ^0, namely, 

2q^ir{R/'-Ro')e^ (68) 



re dr 



and upon reducing this expression to an integrated form, employing a process identical with 
that used for reducing equation (57) , there results as the expression for efFectivemess of a cir- 
cumferential fin, 

^ 2flc i J j (iaR c) fli (iaR't) -i Ji (iaR's) Hi (idRo) /ftQ^ 
{R't'-Ro') Jo iiaRo) H, {mR7)-Jr {iaR/) Ho ^''''^ 

It is convenient to have the result stated in terms of the fin width and ratio of the inner 
and outer fin radii. 
Let 

P-^, (70) 
w' = Rt'~Rc 



HEAT CONDUCTION m THE FINS OF AIB-COOLED ENGINES. 

From these definitions, it follows that 
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1-P 



2 a.(ia.',.^)^,(;-;)_^/wx 



(71) 
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tanhaw' . , % ion 
where = w+ and a - y ^g. 



IV. RECAPITULATION OF MATHEMATICAL DERI VAT, ONS^CONCLUSIONS 

Of heat conduction lead to the expression fundamental mathematics 

"aw;' 
= width of fin. 
^ = thickness of fin. 

h dissipnted. »<i.we.n tb. fin surface .„d the Mr streun into which the halt 

bKi of the fin; i. e., the te^oLTm^rf.J^^^' ',' W»nti«J "ilh that of tl« 

the fin. ' """'"^ »3'>"'l"w«ll along the line of. tt<«am»ntrf 

•ssu-^riot «" «"W >» "SO„-"y ».ot onder the following 

the expend end. i. .n i/applS Z^'^'"' 
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<3, Co„dit.o,» at .he expend edg, .uch '^SZu^ 
by adding to the J-;^^^ tughX% blanketed. 

" VCtL"d1. pi'aLg .t theM of the Bn i. that of nnrfonn, con- 



8tant temperature. 

(5) The fin thickness is uniform. 



(6) The fin is of rectangular contour. setting up the equations 

By removing restrictions expressed m --7^^74(1;",^^ rlmdning'f our con<Utions, a 
to express exactly the th^'^^l^ the co^^^^^^^ to apply to the function/' . 

solution may be obtained which "^'^r*;;*^^^'^^^^^^^ making assumptions (1) and (3). 
quoted above to take account obvious that proper procedure 

Were this exact solution somewhat "^f ^. "^'^^^f^,^ 

would be to employ it directly in computaUons ^^^^ S;'^;^^^^ .nsuited to repeated 

an mexact formula. However it is found to and physical properties 

use. Fortunately, it proves that for -jf ^^^^^^^^ Sb within 1 per Lt, and for the 
likely to charactenze an engine cyl;nder fin the co rection . ^^^^^^^^^ ^^^^^^^^ 

usual present^ay ^1^^^ mathematics to practical problems 

comparison with other errors inherent m api»y|ng . ^ ^ Ucations of these mathe- 

In view of the foregoing it is quite J^f^^^^'^^^fXi™^^^^^ the hypothetical con- 

matical developments to neglect - -^^^^3^^^^ ^^^^^^^^^^^ which do occur. It is 
ditions outlmed m ^ssump ions (1 and^(3) and th^ 

entirely satisfactory to employ the formula basea on in« 1 , ^^^in satis- 

very exactly the actual fin P^^--^-. . ^^^^ (^^^^^^^^^^^ fins 
factory limits for any numerical work ^ th radxul ^^^^^^^^ connection with assumption (6). 
any error due to the assumption merges into X^rnlles ^^^^^^ requiring any modifica- 
The linutation -r^f ^^^^^^^^^^^^ ;lX::te:Utions\utliIed in the first 

tion of the function / . Provided that ^^^^^^ . ^ . (5) ^nd (6), it is a compara- 
three assumptions and assume the rectangulY 00^^^^ effectiVeness when the 

tively simple mathematical problem to derive the expressio 

The assumptions numbered (5) and W are l"^/^"''" ; ' deviation from the 

practice, and it is therefore of P''-^ -P^^^^^,!^^^^^^^ based upon 

Conditions described in these ^J^'^ J^^^^^^^^^ which is straight in iis length 

them becomes absurd. The exact solution lor 1 circumferential fan is 

coordinate (i. e., a trapezoid P™^^^^^^^^^^^^ ^ is the result of either solution a 

not a problem offering serious ^-'^f^^'^'^l^^lf^^^^^^ But the functions 

prohibitively complicated expression for u^^^^^^ ,.^^.^^,1 ,,f,^,liar 

occurring (Bessel functions of both kinds ^^^^^ToJ^h^values for numerical work are not 
to others than specialists - ^^^^^^^^^f^^^^^^^^ therefore, to give in this paper 

always conveniently accessible. It has seemea v«ry , function 

values of the differences between the ^^^^^"^^ ^'^^^^^ 

tanh aw'law', so that the latter expression might always be used as tlie 
tion and corrections appUed for its error expressed in terms of the fin width (3) 

For the trapezoid section prism fin the results express ^ 
and a taper ratio defined as the ratio of the thickn s ^^J''^ J'P/^-^^^;^'^ ^ sharp-edge 
uniform Lper being assumed^ ^^nTas x'- ^ Si to wMt^ i i co^ttted to by adding 
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on the family of curves forming Figure 9 or those forming Figure 10 and read oflF a quantity to 
be added to the approximate effectiveness. This will give the fin effectiveness as if computed 
from the exact Bessel function formula. The Hmitations to accuracy of the method are those 
imposed by the curves of Figures 9 and 10, which have been computed to the highest accuracy 
convenient for the tables of Bessel functions, etc., usually at hand. Generally speaking, this 
was within 1 part in 1,000 for each individual interpolation made, and for the result when all 
factors arid terms are brought together it is quite certain that the corrected/' is reliable within 
1 per cent, probably within a few tenths of 1 per cent. 

For a circumferential fin, the process suggested is very similar to that just outlined for the 
correction on account of taper. Results are expressed in terms of fin width w, which is the dif- 
ference between the outer and inner fin radii and a ratio p of the inner radius of curvature to the 
outer. It is obvious that for very small values of p, namely, the configuration approaching a 
circular plate with no hole (p = 0), it would be absurd to use a formula based upon a long, rec- 
tangular fin and the ^'corrections'' to such formula which have any real validity as corrections 
are limited to the larger values of p. Between the values p = 1 and p = 0.5 the method is appli- 
cable, but toward the lower value of p the corrections become large and likewise less certain. 
After increasing iv by half the fin thickness to u/, compute tanh aw' jaw' and then in terms of 
this value of approximate effectiveness and p, the curvature ratio, make use of Figure 11 or 12 
to ascertain the correction to be added to the approximate value to obtain the true value of fin 
effectiveness. It may be noted that this correction is always negative; i. e., a circumferential 
fin is less effective than the value computed from tanh av//aw\ 

Summing up the foregoing paragraphs with respect to the diflference between conditions 
which would meet assumptions (1) to (6) and the conditions which actually prevail, it is to be 
noted that only the differences concerned in (5) and in (6) have appreciable effect upon com- 
putations for the fins of internal-combustion engines. In case both assumptions are violated 
at once, namely, a tapering fin of circumferential type, two corrections may be applied to the 
approximate function, with a somewhat less degree of accuracy than pertains to either correc- 
tion alone. The correction for taper was determined on the hypothesis of no circumferential 
curvature and the correction for curvature on the hypothesis of uniform thickness, whence it is 
clear that if both factors are concerned the corrections applied by this method are not exact 
The deviation is a second-order error and is usually too small to be significant in industrial 
applications of such a computation. 

V. EXAMPLES OP COMPUTATIONS. 

For the purpose of illustrating the ease of using the formulas whose derivation has been 

discussed above, a few examples are appended. These are selected from two well-known 
aviation engines— the Gnome, which has steel fins, and the Lawrance, with aluminum fins. 
A long, straight fin of the same width and equivalent uniform thickness is used for an initial 
example, followed by approximating to (a) the taper (6) the annular shape of the real engine 
fin. For physical interpretation of the results of the computations, the definition of effective- 
ness must be borne in mind. An effectiveness of 85 per cent means that each small area of fin, 
say, 1 cm.^ or 1 sq. in., is equivalent in heat dissipating power to 0.85 as much area all at a 
temperature identical with that of the engine cylinder wall in that vicinity. 

The computations illustrate in a convincing manner the relative unimportance of high 
thermal conductivity for fin metal. While the difference in effectiveness of the steel aiid alumi- 
num fins is quite appreciable, nevertheless it may be noted that the effectiveness of steel fins, 
very thin and yet reasonably wide, is high enough that for the conditions assumed in these 
examples it is comparable to that of aluminum fins having five times the thermal conductivitv 
of the steel ones. Thus, no great incentive exists, on this score at least, to employ metals of 
very high thermal conductivity. 
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Example 1. — Si-eel fin, long, redangvlarj uniform thickness. 

Width, 1.60 cm. 
Thickness, 0.08 cm. 

Thermal conductivity of steel, assume 0.10 cgs. units. 

Surface dissipation coefficient— assume g = 0.003 cgs. units, corresponding probably to 
free-air speed in the vicinity of the engine of 40 to 50 meters per second, 90 to 110 miles per 
hour. 

t 



a 



w 



t(?+|= 1.60 + 0. 04«1. 64 



^ ^kt V 0.10X0.08 ^ 



aw' -1.420 
tanh aw' - 0. 8896 

tanha^^^ ^26 
aw 

The effectiveness of such a steel fin is thus about 63 per cent. 

Example 2.— Steel fin, long, redangidar, tapering. 

Width 1.60 cm. 

Thickness, 0.05 cm at tip and 0.11 cm at root. 
Thermal conductivity of steel, assume 0.10 cgs units, 
gf, assume 0.003, as in Example 1. 

w'=-w + ^i = l .600 + 0.025-1 .625 cm. 



l2q /" 2X0.003 ^ RfiS 



aw' = 1.408 
tanh aw' = 0.8870 

tanW'^Qg3Q 

aw 

The taper factor X = <i/<m is 0.05/0.08 = 0.625 and from Figure 10 the correction for a taper 
factor 0.63 and approximate effectiveness 0.63 is O.O34. This correction added to 0.630 gives 
0.66,. 

The efiFectiveness of such a steel fin is therefore 66 per cent. 

It is perhaps worth while to illustrate the computation of effectiveness of a tapenng fin 
directly from the exact equation (58) which has furnished the basis of the corrections plotted 

as Figure 10. ij. fa ,) 77, fa) -iJi fa ) (W 

J'' 2qw' ^= 77^ {in,) Jo fa) +iJi fa) Hlo fa) 
Th? half angle a ( = tan a) of the wedge is 0.0187. 
From this, by equations (52) 

Mi = 2.90 
M, = 4.34 
J„ = 15.17 
"i/, (ini) = 3.613 
iJi (W = - 13.30 
iHo (in,) = 0.00487 
(iAi.)= -0.0288 
H, (ipj) = 0.00540 

2^-0-^1 
(-3.613) (-0.00540) -(-13.30) (-0.0288) 
(-0.0288) (15.17) + (-3.613) (0.00487) 
/= 0.665 



" Jaholce »nd Emde, Fnnktioneotaleta, pp. 130-135, 1809 edition. 
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Example Z—Stedjiny circumferential ^ thickness uniform. 

Width, 1.60 cm. 
Thickness, 0.08 cm. 
Inner radius, 5.65 cm. 

Thermal conductivity of steel, 0.10 cgs. units. 
2 = 0.003 cgs. units. 

From Example 1 the first approximation to effectiveness is 0.626. The outer fin radius, 
7.25 cm. plus half the thickness, is 7.29, whence the circumferential curvature factor o-^RJR\ is 
5.65/7.29, or 0.775. 

From Figure 12 the correction for a circumference factor of 0.78 and approximate effective- 
ness 0.63 is -0.032. Adding this to 0.626 gives 0.594. 

The effectiveness of this fin is therefore approximately 59 per cent. 

To illustrate the computation of the above example directly from the exact solution for an 
annular fin, equation (71), instead of through the medium of a correction curve based on this 
equation, the following figures are summarized: 

From example 1, aw' = 1.420 
P = 0.775 
1-p -0.225 

aw'7-^ = 4.89 
1-p 

f<=6.32 

1 -P 

Jo (^mw' J ^^^ = 24.69 
ij(mw'^^^^ -22.00 

\1^/ -84 (estimated) 
ilio ^ia^Y^-^- 0.00265 

( r^)== -0-000605 
//i (^^»^'rz^)= -0.00291 



aS^^lrrpJ^ 0.615 

f-Q6H^~^^'^^^ (-0.00060 5) -( -84) (-0.00291) ^ 
J u.Dio (24.697 ( - 0.000605) T(^:^84r(0.00265) ~ ' 

Note.— The value 84 for ^ (i 6.32) can not be ascertained with any precision, but by using 
an identical value in both numerator and denominator, the accuracy of computation is not 
vitiated more than 1 or 2 parts in 500 by the probable error. The agreement of 0.594 with 0.597 
is well within the limit to be expected in the use of Figure 12. 

With a steel fin having the taper of Example 2 and the curvature of Example 3, with the 
remaining characteristics the same as those taken for all three examples, the effectiveness would 
be approximately 0.594 (Example 3) plus 0.034 for taper, or 0.63. 
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(ft) Computation of optimum dimensions of fins « to meet any specified relations which are 
not mutually inconsistent. An example of this is minimum weight compatible with a given 
cooling requirement, or head resistance, or combination of the two. 

In applying the equations to compute optimum dimensions, it must be borne in mind that 
too much must not be expected in the way of results. There are always a number of conditions 
to be met which arc conflicting and for which a rigid mathematical specification of all would 
render the problem unsolvable. In fin design these include maximum cooling power, minimum 
weight, minimum head resistance to the air stream, adequate strength to withstand crushing 
under rough handling, a choice of metal and geometrical form consistent with the possibiUty of 
manufacturing with reasonable convenience and without prohibitive cost, etc. The engineer 
has no grounds for expecting mathematics to furnish a single inviolable solution for the optimum 
dimensions of the fin that will meet best such an array of specifications, but he does have the 
right to expect mathematics to furnish him definite relations in which he can weigh the various 
factors. Then, according to his judgment of relative importance of such factors, he can select 
the design which he considers best. It is the purpose of this paper to supply definite relations 
respecting the effectiveness of cooling for fins of ordinary type. 
Bureau of Standards, 

WasJiingtonj D, C\, December 1, 1921. 

""^xamples^oT^immQ diniension c^xlculatioa seem best deferred to a separate paper so as to permit of more detailed discussion than 
can be included here. 
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Poaitive directions of axes and angles (forces and moments) are shown by arrows. 





r orce 
( parallel 
to axis) 
symbol. 


Desiignation. 


Sym- 
boL 


Longitudinal 


A' 


X 


Lateral 


r 


Y 




z 


Z 



Moment about axis. 




Angle. 



Designa- 
tion. 



Sym- 
bol. 



roll ! * 

pitch I e 

yaw j ^ 



Velocities. 


Linear 

(oompo- 
neut along 
axis). 


Angular. 


u 

V 

w 


V 
<1 
r 



Absolute coefficients of moment 
^'^qb~S ^--gc5 ^^^ifS 



Diameter; D 

Pitch (a) Aerodynamic pitch, 

(b) Effective pitch, 

(c) Mean geometric pitch, pg 

(d) Virtual pitch, p^ 

(e) Standard pitch, 2>s 
Pitch ratio, pfD 

Inflow velocity, F' 
Shpstroam velocity, F» 



Angle of set of control surface (relative to 
neutral position), 5. (Indicate surface by 
proper subscript.) 

4. PROPELLER SYMBOLS. 

Thrust, T 
Torque, Q 
Power, P 

(If ''coefTicients" are introduced all units 
used must be consistent.) 
Efficiency T V/P 
Revolutions per sec, n; per min., N 

\2'irrn/ 
5. NUMERICAL RELATIONS. 



Effective helix ande ^ = tan" 



1 BP -76.04 kg. m/sec.-550 lb. ft/sec. 
1 kg. m/sec. = 0.01315 IP 
1 nii/hr. = 0.44704 m/sec. 
1 m/sec. = 2.23693 mi/hr. 



1 lb. -0.45359 kg. 

1 kg. -2.20462 lb. 

1 mi.=3lG09.35 m. = 52S0 ft. 

1 m. =3.28083 ft. 
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